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We present applications of variational - wavelet approach to three different models 
of nonlinear beam motions with underlying collective behaviour: Vlasov-Maxwell- 
Poisson systems, envelope dynamics, beam-beam model. We have the represen- 
tation for dynamical variables as a multiresolution (multiscales) expansion via 
high-localized nonlinear eigenmodes in the base of compactly supported wavelet 
bases. Numerical modelling demonstrates formation of coherent structures and 
stable patterns. 



1 Introduction 



In this paper we consider the applications of a new numerical-analytical 
technique which is based on the methods of local nonlinear harmonic anal- 
ysis or wavelet analysis to three nonlinear beam/accelerator physics prob- 
lems which can be characterized by collective type behaviour: some forms of 
Vlasov-Maxwell-Poisson equations [1], RMS envelope dynamics [2], the model 
of beam-beam interactions [3]. Such approach may be useful in all models in 
which it is possible and reasonable to reduce all complicated problems re- 
lated with statistical distributions to the problems described by systems of 
nonlinear ordinary/partial differential equations with or without some (func- 
tional)constraints. Wavelet analysis is a relatively novel set of mathematical 
methods, which gives us the possibility to work with well-localized bases in 
functional spaces and gives the maximum sparse forms for the general type 
of operators (differential, integral, pseudodifferential) in such bases. Our ap- 
proach is based on the variational- wavelet approach from [4]-[14], which allows 
us to consider polynomial and rational type of nonlincaritics. The solution 
has the following multiscalc/multircsolution decomposition via nonlinear high- 
localized eigenmodes 



u(t,x)= a t3 W(x)V 3 (t), 



V\t) = Vp slow (t) + ]T Vffrlt), u\ ~ T 



(1) 
(2) 



i>N 
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'i 



2->. 



(3) 



j>M 



which corresponds to the full multiresolution expansion in all underlying 
time/space scales (x are the generalized space coordinates or phase space 
coordinates, t is time coordinate). Formula (Q) gives us expansion into the 
slow part Mn°m an d f as t oscillating parts for arbitrary N, M. So, we may 
move from coarse scales of resolution to the finest one for obtaining more 
detailed information about our dynamical process. The first terms in the 
RHS of formulae (l)-(3) correspond on the global level of function space de- 
composition to resolution space and the second ones to detail space. In this 
way we give contribution to our full solution from each scale of resolution or 
each time/space scale or from each nonlinear eigenmode (Fig.l). The same 
is correct for the contribution to power spectral density (energy spectrum): 
we can take into account contributions from each level/scale of resolution. 
In all these models numerical modelling demonstrates the appearence of co- 
herent high-localized structures and stable patterns formation. Starting in 
part 2 from Vlasov-Maxwell-Poisson equations, root-mean-square (RMS) en- 
velope dynamics and beam-beam interaction model we consider in part 3 the 
approach based on variational- wavelet formulation. We give explicit represen- 
tation for all dynamical variables in the base of compactly supported wavelets 
or nonlinear eigenmodes. Our solutions are parametrized by solutions of a 
number of reduced algebraical problems one from which is nonlinear with the 
same degree of nonlinearity and the rest are the linear problems which corre- 
spond to particular method of calculation of scalar products of functions from 
wavelet bases and their derivatives. In part 4 we consider numerical modelling 
based on our analytical approach. 

2 Collective models 

2.1 Vlasov-Maxwell-Poisson Equations 

Analysis based on the non-linear Vlasov-Maxwell-Poisson equations leads to 
more clear understanding of the collective effects and nonlinear beam dy- 
namics of high intensity beam propagation in periodic-focusing and uniform- 
focusing transport systems. We consider the following form of equations 
([1],[2] for setup and designation): 
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Figure 1. Multiscalc/cigcnmodc decomposition. 



(J?+!?>=-t J dp » dp « fb > (5) 

J dxdydp x dpyf b = N b (6) 

The corresponding Hamiltonian for transverse single-particle motion is given 
by 

H(x,y,p x ,p y ,s) = i(p2 + p 2 y ) + ^[k x (s)x 2 (7) 

+ky(s)y 2 } +H 1 (x,y,p x ,p y ,s) + ip(x,y,s), 

where Hi is nonlinear (polynomial/rational) part of the full Hamiltonian. In 
case of Vlasov-Maxwell-Poisson system we may transform (4) into invariant 
form 

^ + [f,H} = 0. (8) 

2.2 RMS Equations 

We consider an approach based on the second moments of the distribution 
functions for the calculation of evolution of RMS envelope of a beam. The 
RMS envelope equations are the most useful for analysis of the beam self- 
forces (space-charge) effects and also allow to consider both transverse and 
longitudinal dynamics of space-charge-dominated relativistic high-brightness 
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axisymmetric/asymmetric beams, which under short laser pulse-driven radio- 
frequency photoinjectors have fast transition from nonrelativistic to relativis- 
tic regime [2]. Analysis of halo growth in beams, appeared as result of bunch 
oscillations in the particle-core model, also are based on three-dimensional 
envelope equations [2] . We can consider the different forms of RMS envelope 
equations, which are not more than nonlinear differential equations with ra- 
tional nonlincarities and variable coefficients from the formal point of view. 
Let f(x\,X2) be the distribution function which gives full information about 
noninteracting ensemble of beam particles regarding to trace space or trans- 
verse phase coordinates (xi,X2)- Then we may extract the first nontrivial 
effects of collective dynamics from the second moments 



Q 2 



>= J J XiXjf{xi,Xj)&Xi&Xj (9) 

RMS emittance ellipse is given by s xrms —< x\ >< x 2 - > — < XiXj > 2 (i ^ 
j). Expressions for twiss parameters are also based on the second moments. 
We will consider the following particular cases of RMS envelope equations, 
which describe evolution of the moments (9) ([2] for full designation): for 
asymmetric beams we have the system of two envelope equations of the second 
order for a Xl and a X2 : 

< + <4 ^ + ^ 2 ^ - I/(Io(a Xl + ^ 2 ) 7 3 ) + e 2 nx Jaiy, (10) 



2 



The envelope equation for an axisymmetric beam is a particular case of pre- 
ceding equations. Also we have related Lawson's equation for evolution of 
the rms envelope in the paraxial limit, which governs evolution of cylindrical 
symmetric envelope under external linear focusing channel of strenght K r : 



where K r = —Fr/r^jmc 2 , (3 = Vb/c = \J\ — 7~ 2 . According [2] we have 
the following form for envelope equations in the model of halo formation by 
bunch oscillations: 

Y + kl(s)Y-f^-i =tl , (12) 
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2 + ^-7^-^ = 0, 

where X(s), Y(s), Z(s) are bunch envelopes, £, X ,S, V , £,z = F(X,Y,Z). 

After transformations to Cauchy form we can see that all these equations 
from the formal point of view are not more than ordinary differential equations 
with rational nonlinearities and variable coefficients Also, we may consider 
regimes in which 7, 7' are not fixed functions /const ants but satisfy some 
additional differential constraints/equations, but this case does not change 
our general approach of the next part. 



2.3 Beam-beam modelling 

In A. Chao e.a. model [3] for simulation of beam-beam interaction the ini- 
tial collective description by some sort of equation for distribution function 

f(s,x,p) 

| +P |-( fc ( S ),-F(,, s) )| = (13) 

is reduced to Fockkcr-Planck (FP) equation on the first stage and to very 
nontrivial dynamical system with complex behaviour, 



d,2 1 -*ds 0ia\a\ 

^=T k a k (l-alat), (14) 

which solution gives the parameters of enveloping gaussian anzatz for solution 
of FP equation, on the second stage. From the formal point of view equations 
(14) are particular case of system (12). 

3 Rational Dynamics 

After some anzatzes [15] our problems may be formulated as the systems of 
ordinary differential equations (cases 2.2 and 2.3 (system (14)) above) 

Qi{x)—£j- = Pi(x,t), X=(X!,...,X n ), (15) 

i = l,...,n, max deg Pi = p, ma,xdegQi = q 

i i 

or a set of such systems (cases 2.1, 2.3 (full equation (13)) above) corre- 
sponding to each independent coordinate in phase space. They have the fixed 
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initiator boundary) conditions 2^(0), where Pi,Qi are not more than poly- 
nomial functions of dynamical variables Xj and have arbitrary dependence 
of time. Because of time dilation we can consider only next time interval: 
< t < 1. Let us consider a set of functions 

= Xi-^iQiVi) + PiVi (16) 

and a set of functionals 

Fi{x) = f $i{t)dt - Own In, (17) 
Jo 

where yi(t) (j/i(0) = 0) are dual (variational) variables. It is obvious that the 
initial system and the system 

Fi(x) = (18) 

are equivalent. Of course, we consider such Qi(x) which do not lead to the 
singular problem with Qi(x), when t — or t = 1, i.e. Q,(a;(0)), Qi(x(l)) 
oo. 

Now we consider formal expansions for Xi , y, : 

Xi(t) = Xi (0) +^A^fe(t) yj(t) =^))>rW, (19) 

k r 

where tpk(t) are useful basis functions of some functional space (L 2 ,X P , 
Sobolev, etc) corresponding to concrete problem and because of initial condi- 
tions we need only tpk(0) = 0, r = 1, N, i = 1, n, 

A = {A i } = {A[} = (A^A?,..,Af), (20) 

where the lower index i corresponds to expansion of dynamical variable with 
index i, i.e. Xi and the upper index r corresponds to the numbers of terms 
in the expansion of dynamical variables in the formal series. Then we put 
( |l9| ) into the functional equations (|l^) and as result we have the following 
reduced algebraical system of equations on the set of unknown coefficients A^ 
of expansions (|l^) : 

L{Q ij ,\,a I ) = M(P ij ,\,j3j), (21) 

where operators L and M are algebraization of RHS and LHS of initial problem 
(p"5|), where A ( pp| ) are unknowns of reduced system of algebraical equations 
(RSAE)©. 

Qij are coefficients (with possible time dependence) of LHS of initial sys- 
tem of differential equations ([l5]) and as consequence are coefficients of RSAE. 
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Pij are coefficients (with possible time dependence) of RHS of initial sys- 
tem of differential equations ([l5]) and as consequence are coefficients of RSAE. 
I = iq+2), J = (iij jp+i) are multiindexes, by which are labelled aj 

and (3i are other coefficients of RSAE (^lj) : 

fr = {&Wp+J= / II ( 22 ) 
l<?'fc<p+l 

where p is the degree of polinomial operator P fllq ) 

aj = {a h ...a iq+2 } = cp h ...tp ia ...(p^ +2 , (23) 



where q is the degree of polynomial operator Q (|15|), ie = (1, <? + 2), i/3'j s = 
d<ft 3 /di. 

Now, when we solve RSAE ( pl| ) and determine unknown coefficients from 
formal expansion ( |l9| ) we therefore obtain the solution of our initial problem. 
It should be noted if we consider only truncated expansion (|l9|) with N terms 
then we have from ( pl| ) the system of N x n algebraical equations with degree 
I = max{p, q} and the degree of this algebraical system coincides with degree 
of initial differential system. So, we have the solution of the initial nonlinear 
(rational) problem in the form 

N 

x l (t)=x l (0)+Y / ^Xk(t), (24) 
fc=i 

where coefficients are roots of the corresponding reduced algebraical (poly- 
nomial) problem RSAE ( |2l| ) . Consequently, we have a parametrization of so- 
lution of initial problem by solution of reduced algebraical problem ( pl|). The 
first main problem is a problem of computations of coefficients ai ([2 3D , (3j 
( p2| ) of reduced algebraical system. These problems may be explicitly solved 
in wavelet approach [4]- [6]. The obtained solutions are given in the form 
( p4| ) , where Xk (t) are basis functions and X l k are roots of reduced system of 
equations. In our case X^it) are obtained via multiresolution expansions and 
represented by compactly supported wavelets and are the roots of corre- 
sponding general polynomial system (^T|). Because afhne group of translation 
and dilations is inside the approach, this method resembles the action of a 
microscope. We have contribution to final result from each scale of resolution 
from the whole infinite scale of spaces. More exactly, the closed subspace 
Vj(j £ Z) corresponds to level j of resolution, or to scale j. We consider a 
multiresolution analysis of L' 2 (R n ) (of course, we may consider any different 
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functional space) which is a sequence of increasing closed subspaces Vj : 

...V-2 c V_i c y c Vi c V 2 c ... 

satisfying the following properties: let Wj be the orthonormal complement of 
Vj with respect to Vj + \ : Vj+i = Vj Wj, then 



L 2 (R) = Vo0V^, (25) 

This functional space decomposition corresponds to exact nonlinear eigen 
mode decompositions (2), (3). It should be noted that such representa- 
tions give the best possible localization properties in the corresponding 
(phase) space/time coordinates. In contrast with different approaches formu- 
lae (l)-(3) do not use perturbation technique or linearization procedures and 
represent dynamics via generalized nonlinear localized cigenmodes expansion. 
So, by using wavelet bases with their good (phase)space/time localization 
properties we can construct high-localized coherent structures in spatially-ex- 
tended stochastic systems with collective behaviour. 

4 Modelling 

Resulting multircsolution/multiscale representations for solutions of equa- 
tions from part 2 in the high-localized bases/eigenmodes are demonstrated 
on Fig. 2-Fig. 7. This modelling demonstrates the appcarcncc of stable pat- 
terns formation from high-localized coherent structures. On Fig. 2, Fig. 3 we 
present contribution to the full expansion (l)-(3) from level 1 and level 4 of 
decomposition (25). Figures 4, 5 show the representations for full solutions, 
constructed from the first 6 cigenmodes (6 levels in formula (25)). Figures 6, 
7 show stable patterns formation based on high-localized coherent structures. 
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Figure 2. Eigcnmode of level 1 




Figure 3. Eigcnmode of level 4 
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Figure 4. Appcarcncc of coherent structure 




Figure 5. Six-eigenmodcs decomposition 
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Figure 6. Stable pattern 1 
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